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CHAPTER  I 


THE  REVERBERATION  MODEL 

I.  INTRODUCTION 

In  order  to  predict  the  performance  characteristics  of  an  active 
sonar  in  a  particular  ocean  environment  one  needs  a  certain  amount  of 
knowledge  concerning  the  nature  of  the  background  interference.  Since 
reverberation  is  a  prime  source  of  Interference  and  since  it  is  also  a 
random  process,  the  theoretical  description  must  be  in  terms  of  co- 
variances,  Joint  probability  densities  and  other  such  probabilistic 
quantities.  The  following  sections  develop  one  possible  description 
of  the  reverberation  in  these  statistical  terms. 

Before  beginning  the  detailed  description  of  reverberation  it  is 
advisable  to  consider  briefly  the  assumptions  which  arc  used  in  this 
description.  The  first  is  that  the  reverberation  noise  can  be  ade¬ 
quately  specified  for  engineering  purposes  by  knowledge  of  the  covariance 
function  only.  The  implication  of  this  assumption  is  that  a  knowledge 
of  the  covariance  and  corresponding  power  density  spectrum  is  sufficient 
for  both  engineering  purposes  and  investigations  of  the  statistical 
nature  of  the  reverberation  noise  amplitude.  The  second  assumption  is 
that  the  properties  of  the  various  scattering  media  present  in  the  ocean 
are  adequately  described  by  the  various  experimental  results  reported  in 
the  literature.  Finally,  it  is  assumed  that  the  Fraunhofer  diffraction 
model  used  to  characterize  array  patterns  and  the  ray  tracing  method 
for  calculation  of  the  geometrical  propagation  loss  are  both  valid 
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approximations  for  the  purposes  of  the  reverberation  description.  The 
next  step  Is  to  outline  the  method  by  which  the  covariance  of  the  rever¬ 
beration  process  is  obtained. 


II.  THE  COVARIANCE  FUNCTION 
1.  Preliminary  Remarks 

The  primary  emphasis  In  the  following  sections  will  be  placed  on 
finding  the  covariance  function  resulting  from  scattering  from  a  deep 
scattering  layer.  The  results  thus  obtained  for  volume  scattering  can 
then  be  easily  extended  to  the  case  of  surface  and  bottom  scattering. 

The  method  used  to  obtain  the  covariance  function  may  be  outlined 
as  follows: 

a.  The  transmitted  signal  is  first  written  in  terms  of 

its  time  and  spatial  dependence  utilizing  the  Fraunhofer 
Diffraction  model  for  the  array  and  the  geometrical 
spreading  loss  expression  from  Chapter  II. 

b.  The  interaction  of  the  outbound  pressure  wave  with  the 
individual  scatterers  is  evaluated  utilizing  the  simplest 
scattering  model  which  will  give  results  consistent  with 
the  measured  results. 


The  returning  waves  from  all  scatterers  are  then  summed 
at  the  receiving  array,  after  correction  for  its  direc¬ 
tional  characteristics.  The  result  is  an  equivalent 
noise  voltage  containing  certain  random  parameters. 


Sirtkur  O  littlrJnr. 


I 

I 

I 

I 

I 

*. 

0 

D 

0 


l 


0 

0 

0 

0 


d.  The  next  step  Involves  writing  the  covariance  function 
for  the  reverberation  noise  voltage,  taking  the  required 
statistical  averages,  and  then  simplifying  the  resulting 
integrals . 

e.  The  final  step  is  to  obtain  the  Fourier  transform  of  the 
covariance  function  which  represents  the  instantaneous 
power  density  spectrum  of  the  reverberation  noise. 

During  the  course  of  the  derivation  outlined  above  certain  mathe¬ 
matical  difficulties  do  arise.  In  order  to  obtain  a  useful  model  these 
difficulties  must  be  circumvented  by  judicious  approximations.  The  most 
significant  approx imat ion  is  implicit  in  the  above  outline  and  involves 
the  fact  that  the  reverberation  noise  is  in  reality  a  non-stat ionary 
random  process.  A  non-stat ionary  random  process  is  one  for  which  the 
covariance  function  is  time  dependent.  In  general,  such  processes  are 
analytically  difficult  to  treat;  however,  there  is  one  approximation 
which,  if  justified,  greatly  reduces  the  mathematical  complexity  of  the 
problem.  In  particular,  if  it  can  be  shown  that  the  random  waveform  in 
question,  say  r(t),  can  be  factored  into  two  parts  as 

r (t)  -  g(t)  n(t)  (1) 

in  which  g(t)  is  a  slowly  varying  totally  deterministic  function  and 
n(t)  contains  all  of  the  remaining  random  parameters  then  the  covariance 
of  r(t)  may  be  obtained  as  follows.  The  true  covariance  is  written  as 
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Rrr  (t  r)  -  r (t)  r  (t  +  ▼) 


in  which  Che  horizontal  bar  indicates  a  statistical  average  over  all 
random  parameters  in  r(t).  Using  the  fact  that  r(t)  is  factorlzable 
the  covariance  may  be  written  as 


Rrr  (t,T)  -  g(t)  n(t)  g(t  ♦  i)  n(t  +  f)  .  (3) 

Since  g(t)  contains  no  random  parameters  it  may  be  brought  outside  of 

the  expected  value  operation  and  R  (t,T)  becomes 

nn 


Rrr<t*T>  *  8(0  +  T>  "(t)  n(t  +  f) 


Sow,  if  g(t)  is  a  slowly  changing  function  of  time  with  regard  to  the 
covariance  function  of  n(t)  alone,  i.e.. 


g(t)  -  g(t  +  t)  <  <  nft)'  -  n (t )  n (t  ♦  ) 


then  for  a  given  instant  of  time,  t  ,  the  true  covariance  function  is 

o 

closely  approximated  by 


R  (t  ,f)  ~  !g(t  )]  n(t  )  n (t  + 
nn  o  o  o  o 


For  the  case  of  active  sonar  reverberation  the  approximation  of 
equation  (6)  is  reasonable  provided  that  the  transmitted  signal  pulse 
resolution  length  is  small  compared  to  the  propagation  range  correspond¬ 
ing  to  the  time  t  .  The  above  considerations  are  shown  in  the  form  of 

o 

a  block  diagram  of  the  calculations  in  Figure  1-1. 
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The  following  model  which  is  derived  in  Chapter  II  is  used  to  cal¬ 
culate  the  geometrical  spreading  loss  along  a  ray  path.  The  loss  is 


given  by 


where 


.  r  cos  t 

f-  -  ue  )  -  - 1 

Xo  °  ftx_ 

*  ht  8 
o  z 


is  the  initial  intensity 


sin  f. 


is  a  unit  distance  from  the  source,  usually  1  yd 


is  the  Initial  ray  angle 


is  the  final  ray  angle 


is  the  horizontal  distance 


In  addition,  with  reference  to  Figure  1-2  the  following  results 
for  Incremental  volume  and  Incremental  area  are  noted. 
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The  ocean  scattering  model  used  In  this  derivation  la  made  up 
of  three  principal  scattering  regions  as  shown  in  Figure  1-3.  The  dis¬ 
tribution  of  scatterers  on  the  surface  and  bottom  and  within  the  deep 
scattering  layer  Is  assumed  to  be  uniform  with  the  following  densities: 

■  average  number  of  scatterers  per  unit  bottom  area 

•  average  number  of  scatterers  per  unit  surface  area 

•  average  number  of  scatterers  per  unit  volume  within 

the  deep  scattering  layer. 

The  basic  model  for  the  individual  scatterers  in  each  region 
is  described  by  the  following  assumptions: 


K 


a.  An  individual  scatterer  is  characterised  by  a  complex 

scattering  amplitude,  whlch  contains  both  the 

magnitude  and  phase  information  for  the  scattered 
wave.  The  subscripts  refer  to  the  physical  location 
of  the  scatterer. 

b.  The  expected  value  of  the  complex  scattering  amplitude 
is  sero: 


E  [Z 


0 


(10) 


c . 


The  returning  waveforms  from  two  scatterers  located 
at  different  points  in  the  ocean  are  uncorrelated. 
Thus 


r~  ~  . 
E  [Z44.  Z.  ] 
ijk  lmn 


0 


(11) 
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<1.  Due  to  the  motion  of  scatterers  on  both  the  ocean 
surface  and  within  the  scattering  layer  there  will 
be  a  Doppler  shift  associated  with  each  scatterer. 
This  Doppler  shift  will  be  a  random  variable  and 
its  effect  on  the  returned  waveform  will  be  multi¬ 
plication  by  a  factor  of  the  form 
J2nf<j(i  Jk)t 

These  properties  arc  illustrated  in  Figure  1-4. 


As  for  the  gross  effects  of  the  scatterers  such  as  mean 
scattered  intensity  and  dependence  on  angle  of  incidence,  it  will  be 
assumed  that  these  characteristics  are  suitably  accounted  for  by  the 
following  models: 

a.  For  surface  scattering  the  Chapman  and  Harris^ 
surface  scattering  model  is  used.  This  model 
predicts  that  the  average  scattering  level  is 


10  log  A#  • 
where 

8 


•i 

3.3  8  log 

-  158  [V  f 


-  42.2  log  8  ♦  2.6 

1/,  -0.58 

J] 


and 

-  gracing  angle  in  degrees 

V  -  wind  speed  in  knts 
w 

f  ■  frequency  in  cps 

A  complete  description  of  this  model  may  be  found 


(12) 


in  reference  1. 


(1)  Chapman,  R.  P.  and  Harris,  J.  H. ,  "Surface  Rackscattering  Strength 
Measured  with  Explosive  Sound  Sources."  Defense  Research  Board, 

Naval  Research  Establishment,  Canada,  NRE  Report  62/3.  (UNCLASSIFIED) 
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The  bottom  scattering  model  is  given  by 

10  log  A  ■  -  28  ♦  10  log  sin  (gracing  angle)  (13) 

5 

+  10  log  sin  (scattering  angle) 

The  volume  scattering  model  assumes  no  dependence 
on  Incident  angle  and  is  characterised  by  a  single 
constant.  For  the  purposes  of  calculating  pre¬ 
dicted  reverberation  a  mean  square  scattering 
strength  of  10  7  per  cubic  yard  was  assumed. 
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4.  The  Covariance  Function 

The  transmitted  signal  will  have  both  a  temporal  and  spatial 
dependency.  In  particular,  If  the  velocity  excitation  applied  to  each 
transducer  is 

V(t)  -  Re  [  v (t )  e  C]  (14) 

where  Re.  ]  indicates  the  real  part  of  the  quantity  in  brackets,  then 
the  corresponding  pressure  excitation  as  seen  at  some  point  in  the  far 
field  of  the  transducer  may  be  written  as 

s. 

C  IT  „ 

p  u.  v  v  (t  -  r  > 

N 

x  exp  j  .Bn  -  •  rn>  ]  (15) 

n*l 

The  parameters  in  equation  (25)  are 

■  observation  corrdlnatcs 

-  density  of  seawater 

*  Propagation  velocity 

■  propagation  path  length 

«  arbitrary  phasing  of  the  nth  transducer 
«  wavelength 

■  unit  direction  vector 

■  location  vector  of  the  nth  transducer 

■  a  constant  related  to  the  particular 

transducer  type 

»  number  of  elements  in  the  transducer 
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SUMMATION  OF  CONTRIBUTIONS  FROM  ALL 
ELEMENTS  GIVES  PRESSURE  AS  SEEN  AT  A 
POINT  (*,  6,4)  IN  THE  FAR  FIELD 


FIGURE  1-5  INTERPRETATION  OF  THE  FRAUNHOFER  DIFFRACTION  MODEL 

FOR  THE  SPATIAL  DEPENDENCE  OF  THE  TRANSMITTED  WAVEFORM 
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p  (t,  ®t.  .k)  -  ziJk  v  (t  -  ^  [9i  -  no  (t  -  ).  or 


t  • J2n  Cf .(ijk)  ♦  f  (8<  ,  f .)]  t 
d  s  1  j 


x  Lt0.  e 


The  last  term  In  equation  (18)  contains  two  terms  representing  two  dis¬ 
tinct  types  of  Doppler  shift.  The  first  type,  f.(ijk),  is  the  Doppler 

a 

shift  due  to  the  scatterer  motion  and  is  a  random  variable.  The  second, 

f  (8  ,  0  ),  represents  the  Doppler  due  to  the  motion  of  the  source.  In 
si  J 


general,  f  (0  ,  0  )  is  of  the  form 
si  J 


f  (8  ,  0.)  ■  2  f  —  cos  8  cos  0 
s  i  J  c  c  i 


where  f^  Is  the  carrier  frequency  and  v^  is  the  ship  velocity. 

Considering  equation  (18)  as  a  pressure  excitation  applied  to 
the  receiving  array,  the  output  noise  voltage  as  seen  at  the  receiver 
input  may  be  written  as 


n  (O  -  S  S  5  (t  ‘  “T  >  V9i 


l  J  k 


2s 

A  *  r  >•  *4] 

O  C 


,  ,  ,  -iM(.dJk)  +  f  (0  ,  t  )]t  <2( 

«  0,  t»,.  «,]  U»j]  .  d  •  1  1 

where  Q  8  ,  0  represents  the  receiving  array  characteristics  and 
R  i  j 

the  result  has  been  sunssed  over  all  possible  scattering  bodies. 

Having  obtained  an  expression  for  the  noise  voltage  the  co- 
variance  function  of  n(t)  is  obtained  as  follows.  Equation  (20)  is 
entered  into  the  defining  equation  for  the  covariance  function  to  give 
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Rnn  (t,  t  ♦  T)  •  I  ^(t)  n  *(t  + 

■  8  (§  z'lJk  v  (t  -  ^  )  5,  C«t  -  no  <t  -  ^  ).  0,  i 

-J2"[f  (ijk)  +  f  (e.  ,  f>.)]  t 

[•,.  »,]  L  t»j]  .  11  ■  1  1 

x  ^  Z  v  (t  +  T  -  — “  )  Q_  .8  -n  (t+T  -  — a  ),  t> 
i  inn  c  i  i  o  c 

i  mn  J 

.  J2nCf  (i ,m,n)  +  f  (8  ,0  )]  (t  +  t)}  (21] 

x  Q  *  [•  .0  ]  L  [0  ]  e  d  8  '  m  ' 

R  i  m  m 

where  the  symbols  ^  ^  indicate  an  expected  value  operation. 

The  problem  now  is  to  accomplish  a  suitable  simplification  of 
equation  (21).  The  first  step  in  this  simplification  is  to  note  that 
by  the  assumption  of  zero  correlation  between  different  scatterers  the 
cross  product  terms  in  equation  (21)  must  be  zero.  Thus,  since 


(*iJk-  C>  ■ 0 


then 

Rnii  (t.  t  .  »>  -  J  E  E  |?i)k|  ~  <t  .  A  ,  ;  *  (t  +  ,  .  A 

« «T  [»1  - n.  <*  -  c  )•  V  -  no  <*  *  ’  -  -T  >• ») 

,  J2nf  (ijk)  t 
x  Dr  [8t,  0^]  L  [Dj ]2  e  d 

J2nf#(8it  Cj)  y  (22) 
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The  next  step  is  to  simplify  the  expression  further  by  evaluating  the 
required  expected  values.  Now  the  only  random  parameters  in  equation 
(22)  are  the  scattering  amplitude  and  the  scatterer  doppler.  Starting 
first  with  the  scatterer  doppler,  the  expected  value  of  this  quantity 
is  written  as 


E  [< 


)2"  f,  (i,  ),  k)  t 


] 


-  2" 


OD 

<• 


J2"f.(ijk)  t 
a 


p(2Hfd)dfd 


(23) 


where  p(2T7f.)  is  the  probability  density  function  for  the  Doppler  shift 
d 

of  the  scatterers.  It  Is  assumed  to  be  the  same  for  all  scatterers  with¬ 
in  a  particular  scattering  layer  or  on  a  particular  scattering  surface. 
The  evaluation  of  equation  (23)  is  particularly  easy  since  its  form 
corresponds  directly  to  the  form  of  the  characteristic  function  of  a 
probability  density  function.  Thus 


E  [< 


)2"f  .(ijk)  t 


]  -  M(t) 


(24) 


where  M(T)  is  the  characteristic  function  of  pif  ,(ijk) j .  From  various 

U 

experimental  measurement s ^  the  form  of  p.f,(l)k).  has  been  found  to  be 

d 

Gaussian  to  a  reasonable  approximation  and  hence 


M(t) 


’2 

a 


(25) 


where  0  is  the  variance  of  p.fJ(l)k)j.  The  range  of  is  between  one 
d  d  d 

and  ten  cycles  per  second  for  carrier  frequencies  from  3.5  to  5  kc  and 

2 

there  is  evidence  that  for  surface  reverberation  o  is  a  function  of  the 

d 

angle  of  incidence  of  the  excitation  wave. 
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The  second  expected  value  operation  to  be  performed  is  that 


pertinent  to  finding  the  mean  square  scattering  strength.  The  required 
integral  is  written  as 


(26) 


where  p  (Z  )  represents  the  probability  density  function  for  the  com- 

1  J  K 

pl ex  scattering  amplitude.  Since  the  scatterers  are  assumed  to  be 
localized  to  certain  regions  of  the  ocean  the  integral  may  be  evaluated 
with  the  aid  of  Bayes'  rule.  Thus 


Pz(Zi  )k)  -  Pz(2,  ^ lA  scatterer)  P.  Scattered 

+  p^fZj^lNo  scatterer)  P.No  scatterer.  (27) 

If  no  scatterers  are  present,  the  scattering  amplitude  is  zero  with 
probability  one  and  therefore  the  second  term  in  equation  (27)  con¬ 
tributes  nothing  to  the  expected  value  Integral.  In  the  regions  in 
which  scatterers  exlxt  the  probability  of  finding  a  scatterer  is 


P  Scatterer]  -  KiV 

Using  these  results  the  mean  square  scattering  strength  is 


(28) 


Uijkl  -KAV 


Jk>  d2l)V 


<Vv 


(29) 


where  By  is  the  mean  square  scattering  strength  per  unit  volume. 
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Inserting  the  results  of  equations  (24)  and  (29)  Into  equa¬ 


tion  (22)  and  using  the  expression 


AV  -  R 


OR 

d0 


sin  0  A8  A0  As 
o  o 


for  the  volume  element  gives 

R  At,  t  +  t)  -  M(T)  21  2  53  v  (t  -  —  )  v*  (t  ♦  t 
nn  2  1  )  k  C 


2s, 


)2"f  (e  ,0  )t 

x  Dr  ,  0  ^ .  L  .  0  j .  e  9  ^ 


2s. 


x  R 


sin  0  A0  A0  As 
00  I  o  o 

o  z 


Carrying  the  sunsnatlons  over  into  integrals  gives 

.2" 


KB., 


Rnn (t  ’  *  +  ’)  "  -p  *<T>  I  d® 


J  d«  J  d0  J  V  *(t  +  ~)  v  *i 


X  V 9  '  °o(t  '  c  K  t]  V  [®  *  n0  (t  + 


?).  *3 


2  )2rf  (9.0)t 
x  Dr  .  0  L  0  e 

I  I 

x  R  —  sin  0,  ds 
of)  1 

0  2 


Further  simplification  is  obtained  by  noting  that 

Ll0]  R 


OR 

*0 


sin  0.  ■  r  cos  0 
1  o  o 


Hence,  equation  (31)  may  be  written  as 
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R  (t.t  ♦  T) 

nn 


-  £>  7*(t  *T 


x  Qj  [•  -  wo(t  -  ^  Qj*  [•  -  fiQ(t  ♦  t  -  ).  t) 

x  D  .8,  0]  L  .0]  r2  cos  0  *g(**0)  d,  (33) 

R  OO 


For  the  general  case  of  time  limited  algnala  the  integrals  of 
equation  (33)  may  be  relatively  easily  evaluated.  The  factors  which 
make  this  determination  possible  are  first  the  spatial  limitations  on 
the  distribution  of  scatterers,  and  second  the  spatial  limitation  on  the 
ensonlflcd  regions  at  any  given  time  due  to  the  time  limited  nature  of 
the  signal.  Both  of  these  properties  will  he  used  to  evaluate  the 
integrals  over  0  and  s. 

The  physical  situation  with  regard  to  the  scattering  layer  is 

shown  in  Fig.  1*6.  At  anv  given  instant  of  time,  say  t  ,  the  ensonlfled 

o 

volume  is  limited  by  the  upper  and  lower  boundaries  of  the  scattering 
layer  and  by  the  pulse  resolution  length.  Since  the  individual  scatterers 
are  uncorrelated  the  integration  over  s  must  be  zero  for  all  s  outside 
the  ensonified  region.  As  a  consequence  of  this  fact  the  integration 

over  0  at  the  time,  t  .  is  limited  to  an  integration  over  0®  as 

o 

defined  by  the  boundaries  of  the  ensonified  volume. 

To  utilise  the  above  results  the  integration  on  s  is  written 

as 
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v  <lo  ■  “T  )  v  UQ  ♦  T  *  “■  )  Qt*  ■  fl©  <‘o  '  ^  >•  •  ' 


<L*  *  -  ..  (t  +  ’  •  ^  iii 


j  *  :‘o  ’  ¥  •  ••  °  *  * 


t  ♦  T  d* 

o  c 


where 


'o  ■  T  •  *•  *  "  *  (to  •  ?  )  ■*  *  o  (to  •  f  >*  * 


Next  a  change  of  variable  is  made  such  that  the  integral  become# 


Hit  -  a,  #.  0  8  *  it  ♦  ▼  -  a,  9,  t  do 

2)o  o 


£  R  (’.«.•) 

2  gg 


where  R  (’.H.O)  la  the  correlation  function  of  the  transmitted  signal 
gg 

Including  anv  effects  due  to  rotating  directional  transmission. 

For  the  case  of  fixed  beam  transmission  the  form  of  the 
complex  dlrectivitv  function  for  the  transmitter  Q^.  <*,#.  is  not  de¬ 

pendent  on  time.  In  addition,  the  excitation  signal  v  (t)  depends 
only  on  time.  Thus,  the  quantity  R^(T.«*.0)  can  be  written  in  a  fac¬ 


torized  form  as 


R  (▼.  #.  t)  -  0  (T)  D  [•.*] 

gg  w  T 
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Using  this  result  in  equation  (33)  and  assuming  for  the  time  being 
that  rotating  directional  transmission  is  not  used  gives 


R  (t,  t  ♦  T)  ■  A  M(f)  t  (f)  r  2  cos  t 
nn  v  vv  o  o 

r  r  • 

y  I  DR(«,f)  D,(«,0)  exp  <J2"f#(e.0)T)d*  (37) 

t  o 

where  we  have  replaced  the  KB^  by  an  experimentally  measured  quantity  A^ 
which  we  will  call  the  scattering  strength. 

The  most  Important  terms  in  equation  (37)  from  the  standpoint 
of  sonar  design  are  the  quantities  H(T)  ,  th*  theta  integra¬ 

tion.  As  noted  above  the  function  M(T)  is  simply  the  Fourier  transform 
or  characterist ic  function  of  the  probability  density  function  describ¬ 
ing  the  random  motion  of  the  scattcrers.  In  addition  to  the  scatterer 
motion  the  function  M(*)  can  also  be  made  to  account  for  random  motions 
of  the  transmitter  and  receiver.  The  overall  effect  of  M(T)  is  to  intro¬ 
duce  Doppler  broadening  to  the  reverberation  return. 

The  correlation  function  of  the  transmitted  signal,  (*VV(T)» 
contains  all  of  the  information  pertaining  to  the  total  energy  and  spec¬ 
tral  content  of  the  transmitted  signal.  For  the  purposes  of  this  dis¬ 
cussion  the  dimensions  of  t>  (T)  will  be  watts  per  square  centimeter. 

w 

The  final  term  in  equation  (37)  involves  an  integration  over 
the  azimuthal  angle,  8,  of  the  receiver  and  transmitter  directivity 
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functions  and  the  exponential  function  containing  the  ships  own  Doppler 
expression.  In  general,  this  integral  must  be  solved  numerically  be¬ 
cause  of  the  complicated  form  of  the  directivity  function.  The  only 
exception  to  this  is  the  case  of  omnidirectional  transmission  for  which 
the  directivity  functions  are  unity  for  all  0  and  0,  in  which  case  the 
integral  reduces  to  a  Bessel  function  multiplying  the  functions  MO) 

and  t  (*)•  The  significance  of  this  integral  lies  in  the  fact  that 
w 

the  ships  own  Doppler  relation  taken  in  conjunction  with  the  directivity 
functions  contributes  to  the  Doppler  spread  of  the  reverberation. 

The  above  considerations  lead  to  the  conclusion  that  the 
dimensions  of  the  covariance  function  must  be  watts  per  square  centi¬ 
meters  when  the  time  shift  T  is  identically  zero.  Thus  the  total  rever¬ 
beration  power  being  returned  at  a  time,  t^,  is  Just  the  evaluation  of 
equation  07)  with  T  set  to  zero.  The  final  step  in  the  analysis  is  to 
obtain  the  spectral  content  of  the  reverberation  by  taking  the  Fourier 
transform  of  equation  (37). 

III.  THE  REVERBERATION  SPECTRUM 

The  Fourier  transform  of  R  (t  O)  is  defined  as 

nn  o 

*  (t  .  f)  -  I  *  <t  ,T)  a‘J2T,fT  dT  08) 

nn  o  J  nn  o 

where  *  (t  ,f)  is  interpreted  as  representing  the  distribution  of  the 
nn  o 

reverberation  power  as  a  function  of  frequency,  f,  at  the  time  t^. 
Inserting  equation  (37)  for  the  covariance  function  gives 
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*  (t  .0  -  \V  2  co»  0  ~  l[0]  dp  DT<0.P)  D10.P) 

nn  o  vo  o  2  I  IT  R 


/'« 


-J2"[f  -  f  (e,P)b 

T)  p  (T)  C  “  dT 

w 


The  final  integral  over  T  is  simply  the  Fourier  transform  of  the  product 
of  two  time  functions  which  is  simply  the  convolution  of  their  respective 


Fourier  transforms.  Hence, 


I  (t  ,f)  •  A  r  cos  0 
nn  o  v  o 


f  r 

I  L.p]  dp  I  Dt(9.P)  Dr(«.P) 


/'-■ 


(a)  t  (f  -  f  (0.0)  -  o)  d- 

vv  s 


where  K  (")  is  the  Fourier  transform  of  M(')  and  C  is  a  dummy  variable 
M 

of  integration. 

So  far  the  discussion  has  been  limited  to  scattering  from  a 
deep  scattering  layer;  however,  it  should  be  pointed  out  that  scatter¬ 
ing  from  the  ocean  surface  and  the  ocean  bottom  may  be  characterized 
by  precisely  the  same  techniques  as  were  used  to  characterize  the  vol¬ 
ume  scattering  problem.  The  answers  are  in  fact  identical  in  form  and 
differ  only  In  terms  of  the  scattering  coefficient  and  the  acatterer 
dispersion  functions.  In  addition,  it  is  noted  that  the  Fourier  trans¬ 
form  of  the  transmitted  signal  autocorrelation  function  may  be  written 


I  CfJ  -IT  A [f] 

W  op 
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where  l  Is  the  peak  transmitted  Intensity  and  t  the  pulse  length, 
o  p 

Thus,  the  product  1  T  contains  the  Information  as  to  the  transmitted 

op 

energy  and  Aif.  gives  the  appropriate  spectral  distribution.  In  the 
following  section  a  complete  sunsnary  of  the  reverberation  equations  is 
given,  including  the  substitution  indicated  by  equation  (<tl). 
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2>  % 

»  \ 

4>  a» 

5)  A, 

6)  c 


Peak  transmitted  intensity  in  watts  per 
square  centimeter 

Pulse  length  in  seconds 

Volume  scattering  strength  per  cubic  yard 
Bottom  scattering  strength  per  square  yard 
Surface  scattering  strength  per  square  yard 
Velocity  of  propagation 


c 

0 


Unit  distance  of  one  yard 


8)  tx 


Angle  of  incidence  of  wavefront  with  scatter* 
ing  surface  or  volume 


9)  fs«».0) 

10)  to 

ll>  ‘o 


Ship's  own  Doppler  function 

Initial  angle  at  the  transmitter  for  the  ray 
which  strikes  a  scattering  boundary  at  time  ft^ 

Round  trip  travel  time  for  the  reverberation 
components 


0 

B 

0 

B 

I 

I 


12)  10] 


Geometrical-spreading  component  of  Intensity 
loss  (one  way) 


13)  f  ;  Frequency  in  c.p.s. 

14)  D^lO.0  and  DR.0,0  ;  Transmitting  and  receiving  directivity 

functions,  respectively 


15)  M(*)  ;  Scatterer  dispersion  function  relating  to 

random  motion  of  the  scatterers 


16)  F  (~)  ;  The  Fourier  transform  of  M(t) 

n 

17)  t  (T)  ;  Autocorrelation  function  of  the  transmitted 

vv 

waveform 


18)  A  [f] 


Spectral  distribution  function  for  the  energy 
in  the  transmitted  waveform 
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b.  Reverberation  Equations 


The  covariance  function  Is  given  by 


R  (t  ,f)  -  A  r2  cos  0  4  M  (t)  0  (t) 

nv  o  vo  o2v  vv 


[  [  J2"f  (e.C)T 

x  I  L  [f]  dt>  J  Dt(O,0)  Dr(9,0)  e 


d9  (62) 


and  the  power  density  spectrum, 

ct 

*  (t  ,f)  •  A  r2  cos  t  I 
nv  o  vo  oo2 


I  L[0] 

"  A 


d0 


r  r 

I  Dt(9.0)  Dr(9,0)  d9  /  Fm<c>*vv 


(f  -  f  (9,0)  -  C)  dc  (63) 


The  only  significant  differences  betweeen  the  various  types  of 
reverberation,  surface  volume  or  bottom,  aside  from  their  respective 

scattering  models  Is  the  effect  form  of  the  characteristic  function  F  (c) . 

M 

For  the  purposes  of  general  systems  analysis  application,  the  assumption 

of  a  gaussian  form  for  F  (o)  Is  reasonable  in  view  of  the  existing 

rl 

experimental  data. 


1-29 


Arthur  D.ltnlrJnr. 


V.  THE  EVALUATION  OF  THE  REVERBERATION  EQUATIONS 

The  numerical  evaluation  of  the  reverberation  equation  begins 
with  an  algorithm  named  T1MBR  whose  purpose  is  to  characterize  a  par¬ 
ticular  ocean  environment  in  terms  of  the  initial  ray  angles  required 
to  propagate  a  ray  to  a  predetermined  point  in  both  space  and  time. 

The  details  of  this  algorithm  are  given  in  Chapter  IV.  The  basic  method 
used  by  the  TIMBR  algorithm  is  essentially  an  exhaustive  characterization 
in  which  rays  are  traced  from  a  source  over  a  very  wide  aperture  and  out 
to  propagation  times  equal  to  or  greater  than  the  maximum  reverberation 
time  of  interest.  A  reasonably  fine  grid  of  initial  angles  is  used  to 
cover  the  aperture  and  the  ray  parameters  such  as  range  and  depth  are 
recorded  at  regular  t  inn-  intervals  as  the  rays  are  traced.  This  results 
in  the  creation  of  a  number  of  TIMBR  arrays  and  corresponding  to  a  par¬ 
ticular  time  and  containing  all  of  the  ray  parameters  for  that  time  for 
all  of  the  rays  in  the  aperture.  Thus  timer  provides  the  ray  parameters 
such  as  range  and  depth  both  as  tabulated  functions  of  initial  angle  and 
as  tabulated  functions  of  time.  Having  obtained  the  TIMBR  description 
of  the  particular  ocean  environment  it  is  then  possible  to  Interpolate 
in  the  TIMBR  tables  by  the  method  described  in  Chapter  IV  to  obtain  the 
integration  ranges  over  the  variable  #>  to  be  used  in  evaluating  equation  (43). 
In  addition,  the  initial  angles  to  be  used  in  the  evaluation  of  the  propa¬ 
gation  loss  are  also  made  available. 

Given  that  the  initial  ray  angles  and  depression  angle  t>  inte¬ 
gration  range  have  been  determined,  the  next  step  is  the  evaluation  of 
the  innermost  signal  convolution  Integral  which  is  accomplished  by  means 
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of  the  fast  Fourier  transform  as  described  in  Chapter  V.  The  next 
integral  is  in  azimuth  over  the  transducer  and  receiver  arrays  and  is 
by  far  the  most  time-consuming  since  the  array  directivity  must  be 
calculated  by  the  Fraunhofer  diffraction  method  once  for  each  step  of 
the  integral.  Finally,  the  integration  over  depression  angle  0  is 
performed  to  give  the  reverberation  spectrum  at  a  given  time. 

The  remaining  chapters  give  complete  details  of  the  methods 
of  performing  each  of  the  calculations  listed  above. 
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CHAPTER  II 


GEOMETRICAL  SPREADING  LOSS  -  GENERAL 


Tlie  spreading  loss  along  a  ray  path  may  be  obtained  by  con¬ 
sidering  the  variation  in  energy  density  in  a  ray  bundle  centered  on 
the  rav  path.  There  are  two  ways  by  which  this  may  be  accomplished. 
First,  consider  the  situation  shown  in  Fig.  II-l  for  which  the  power 
leaving  the  source  through  area  A^  is  given  by 

P  -  I  A  (1) 

o  o  o 


where  l  is  the  initial  intensity  at  a  distance,  r  ,  from  the  source. 

o  o 

At  some  distant  point,  x  ,  Zj,  the  intensity  is  Ij  over  an  area  Aj  such 
that 


P 

o 


I  A 


o  o 


Vi 


(2) 


which  leads  to  the  spreading  loss 


From  Fig.  II-l  note  that  A  is  given  by 

o 

2 

A  ■  r  cos  0  AO  AO 

o  o  o  o 


(3) 


(4) 


where  0  is  the  initial  depression  angle.  The  area  A  at  x  ,  z.  is 
O  III 

given  by 

Aj  •  hw  ■  XjAO  Ax  sinOj  (5) 
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Thus  Che  ratio  of  A^  Co  In  Equation  3  is 


A  r  cos0  A0  A0 

o  o _ o  o 


x,  .'U  ix  slnf^ 


which  in  the  limit  as  ^0  —  0  becomes 

o 


r  “  cos® 


Equation  7  is  perfectly  general  in  that  no  assumptions  have 
been  made  concerning  the  velocity  profile.  However,  if  the  velocity 
profile  is  linear,  then  the  rays  are  known  to  be  arcs  of  circles  and 
the  partial  derivative  Jox/of  j  can  be  evaluated  as  follows.  For  the 
straight  line  velocity  profile  with  gradient  g,  the  horizontal  distance  x 


is  given  by 


cosfM 


s  ■  -  I"  •  —s&  dt  •  -  —  si 
J *  RP  gP 


— *  (sinO  -  sin©  ) 

RP  1  O 


where  p  “ 


and  is  constant  for  a  ray  in  a  constant 


velocity  gradient,  as  a  consequence  of  Snell's  law. 
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Differentiating  Equation  8  gives 


-  (sin*  -  sin*  )] 
3*  3*  ‘g  cos*  1  vo 

o  o  o 


(9) 


■  —j - 5 —  l  (cos*  — - cos*  )  g  cos*  +  (sin*.  -  sin*  )g  sin* 

Z  i.  O  *  O  O  lot 

g  COS  *  o 

o 


The  term  -;/3*  is  obtained  from 
o 


*  *  cos 


.  c ( z )  cos* 
•l  o 


ii- . 

»♦. 


sin; 


£lO 

c 


/  *  sin* 
c(z)  _ _o 


Vl  - 


c(z)‘  COS  * 


c  sin* 
o 


but 


£l£>  . 


cos; 


from  which 


..  cos*  sin* 
li_  .  o 


3*  cos*  sin* 
o  o 


(10) 
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Using  this  result  in  Equation  9  gives 


o  8  cos“»  sin; 

z  °  o 


[sin:  -  sln«) 
o 


cost  sin# 

o 


Using  this  result  In  Equation  7  gives  for  the  spreading  loss 


0 

i  c 
0 
u 


,  2  2 

*1  ro  c°*  >o 

L  "  l - 2 - 

o  x 


This  expression  is  valid  until  the  rav  reaches  a  boundary  between  zones  of 


different  velocity  gradient. 


The  effect  of  the  change  In  velocity  gradient  at  a  zone  boundary  is 
accounted  for  in  the  following  manner.  With  reference  to  Fig.  11-2, 
Equation  7  becomes 


t  r  cost 

L  -  —  -  - - - 2 - 2 _ 

°  x  7J-  (xi  ♦  <V*,>  ♦  (x-x.))  sins 
o 


In  general  for  anv  point  in  the  (n*l)th  velocity  profile  *one. 


o  n+1 
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Differentiating  with  respect  to  # 


J#  <x"*n>  "  2 

°  co>  *o  *n+l 


- — [coMo(co“;n  ir  • COM  U;) 


♦  sin*  (sin*  -  sin#) 
o  n 


] 


From  Snell's  law: 


(15) 


(cos# ) 


Sill  i_ 

c 

o 


o 


c°*» 

o 


-  sin# 


sin: 

o 


sine 

3#  m  cos# _ o 

>*  ’  cos#  sin# 

o  o 


Using  this  result  in  liquation  15  gives 


TT  (x-xn> 
o 


CO*  JoVl 


cos  #  sin# 
n  o 

sin# 


cos  )  sin# 

_ c 

sin# 


♦  sin#  sin# 
n  o 

■  in#  . 

_ o _  _  J _ 

cos#  sin#  *  ini  p  g 

o  n  n*  i 


-  sin# 

(sin# 

n 


sin# 

"J 

-  sin#) 


(x-x  )  sin# 
n _ o 

cos#  sin#  sin# 
o  n 


(16) 
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FIGURE  11-2  PERTAINING  TO  THE  CALCULATION  OF  REFRACTION  LOSS 
IN  A  MULTI  LAYERED  MEOIUM 


Thus  for  the  three  layer  case  given  In  Equation  13, 


L  -  4-  - 


2  2 
r  cos  * 
o  o 


x  i 


sinjj 


(x,-x.)  sin* 

2  1 _ o_ 

•In* ,  slnij 


(x-x  ,)  sin* 
sin*! 


“"ol 

sin* 

J 


(17) 


sin* 


More  generally  for  a  ray  traversing  N  zones  with  a  terminal  point 
In  the  N+l1*'  zone 


L  -  t-  - 


2  2 
r  cos  * 
o  o 


/  s  (x„“X  i*  \  (*"*«.) 

I  - 1  ♦  - i 

\  °  fy  ain:u  ,lnVl  /  altt*S  aint 


(18) 


sin*  ] 


Finally,  if  a  ray  is  specularlv  reflected  from  a  plane  surface,  the 
reflection  must  be  accounted  for  Just  as  though  a  new  zone  had  been 
entered.  Fig.  11*3  Illustrates  this  case  for  surface  reflection. 


The  refraction  loss  in  the  specular  reflection  case  is  calculated 
using 


(x 


n-1 


x  ,)  ♦  (x  -x  .)♦ 
n-2  n  n-1 


(x 


n+l 


Xn)  + 


Just  as  though  the  ray  were  traversing  multiple  layers. 
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CHAPTER  III 


DETAILS  OK  PROPAGATION  LOSS  CALCULATION 


Total  propagation  loss  for  the  selected  rays  consists  of  the 
spreading  loss  defined  by  equation  (18)  in  Chapter  II,  plus  attenuation, 
surface  and  bottom  losses.  Automation  of  the  loss  calculations  requires 
consideration  of  various  special  cases,  such  as  the  tracing  of  vertexlng 
rays.  This  chapter  describes  the  methods  used  to  deal  with  these  situa¬ 
tions.  Details  of  the  computer  programs  are  contained  In  Volume  II  of 
this  report. 


1 .  Ray  Tracing 

The  subroutine  RAYTRK  Is  used  in  both  the  TIMER  and  propa¬ 
gation  loss  calculations.  It  follows  a  ray  leaving  a  starting  point  In 
some  direction,  until  a  specified  amount  of  time  has  elapsed.  TIMER  allows 
the  determination  of  the  angles  at  which  rays  must  leave  the  transmitter 
in  order  to  terminate  on  surface,  bottom  or  volume  scatterers  at  a 
particular  time.  PLOSS  repeats  the  ray  trace  for  each  of  the  selected 
initial  angles  and,  in  addition,  develops  all  of  the  quantities  necessary 
to  evaluate  equation  (16)  In  Chapter  II.  The  ray  tracing  Is  done  In 
steps  along  the  arc  length,  with  the  step  size  adjustable  as  critical 
points,  such  as  vertices,  zone  boundaries  and  end  times  are  approached. 
Newton's  method  Is  used  to  adjust  the  step  size  and  direction  In  order  to 
force  the  rav  trace  to  converge  on  a  specific  point.  The  precision  of 
convergence  Is  a  programmer  controlled  Input. 
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cos; 

c(z) 


cos; 


(2) 


where  cq  Is  the  velocity  at  the  source  and  ;  Is  the  Initial  angle. 

V,  ^  - _ _ t.  .  j 


Thus  the  curvature  may  be  written  as 


ii 


ds  8kPo  "  R^ 


(3) 


where  R^  is  the  radius  of  curvature  in  the  kth  zone.  Hereafter,  the 
subscript  k  indicating  the  kth  zone  will  be  omitted,  and  it  will  be  under¬ 
stood  that  all  equations  hold  within  one  zone. 

The  horizontal  distance  travelled  by  the  ray  is  given  by  * 

*2  “  xi  *  |  cos;  ds 


(4) 


<ing  equation  (3)  to  effect  a  change  of  variable,  we  obtain 


'  V*5  f 


cos;  d*  •  -  ~  (sin;,  -  sin;.) 

RP  2  1 


1 


-  R  (sin;,  -  sin^j) 


(5) 


The  vertical  distance  is  given  by 

.  f*. 


2  2  •  1 


12  1  .? 

j  -  j  sin;  ds  -  -  —  sin;  d; 


*1 


1 


—  (cos;2  -  cosjj) 


-  R  (cos;^  -  cos;j) 


(6) 


*  in  this  and  similar  integrations,  the  integration  is  split  into  two 
rt*  whenever  a  ray  vertexes. 


II 1-3 


3rthur  9Atttlr.)nr. 


Since  the  depth  variable  is  used  by  the  program  as  an  independent  variable, 
equation  (6)  provides  a  means  of  obtaining  the  final  angle  as  follows: 


*2  “ 


cos 


l( 


cos; . 


(7) 


where  p  is  either  ♦  1  and  is  used  to  assign  the  appropriate  sign  to  the 
final  angle  The  method  of  choosing  p  will  be  discussed  when  we  consider 

the  detailed  application  of  the  above  equations. 

The  arc  length  s  is  given  by 


s 


2 


dz 

sin; 


1*2  dz 
*1 


(8) 


Using  Snell’s  Law  and  equation  (3),  note  that 

2  2  2 

cos  $  •  c(z)  p 

o 


(9) 


and 


dz  »  - —  d(c(z)p  ) 
° 


Inserting  these  results  in  equation  (8)  gives 

fZ,  d(c(z)p  ) 


*2  '  S 


.  L_  f*2 

1  KPo  yT  l-c(z)‘p  " 


1 


-  R  (sin  1  c ( z , ) p  -  sin  1  c(z.)p  )  •  0 
i  o  1  o 


(10) 
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The  detailed  application  of  these  equations  requires  consideration  of 
four  possible  cases  as  illustrated  below.  For  the  two  cases  involving  ver- 
texing  rays,  the  calculations  are  done  from  the  initial  point  to  the  vertex 
point  with  complete  evaluation  of  all  integrals.  The  ray  is  then  traced 
from  the  vertex  point  to  the  next  boundary.  This  procedure  is  required 
since  the  vertex  point  is  a  singular  point  and  any  attempt  to  simply 
integrate  through  this  point  will  give  erroneous  results. 


Case  la 


Gradient  <  o  Initial  angle  <  o 

Radius  of  curvature:  R  -  - — 

8P„ 

R  »  o  ° 

Horizontal  distance:  : sine  j i  >  Isini^l 
>  o 

Final  angle:  since  R  >  o  and 
-  Zj  <  o  then 

0  -  -1 


Arc  length: 


sin  1  (c  p  )  <  sin”1  (c,p  ) 
l  o  z  o 


.’.  0  •  -l  in  order  that 


>  o 


NOTE:  Had  the  vertex  depth  fallen  in  the  above  zone,  it  would  have  replaced 
z,  as  the  upper  boundary. 


1 1 1-5 


Arthur  9. little.  Inr. 


Case  lb 


NOTE: 


Gradient  <  o  Initial  angle  >  o 

Radius  of  curvature:  R  ■  -  - —  >  o 

8Po 

Horizontal  distance: 
sin$j  <  sin;^ 

x.  -  Xj  >  o 

Final  angle:  since  R  >  o  and 

Z2  ~  Z1 

z7  -  2^  >  o,  then - - -  <  o 

P  -  +1 

Arc  length:  sin  *  c.p  >  sin  *  c,n  P  -  +1 

1  o  2  o 

such  that  s0  -  Sj  >  o 
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Case  11a  Gradient  >  o  Initial  angle  >  o 

Radius  of  curvature:  R  ■  -  - —  <  o 

8Po 

Horizontal  distance: 

sinij  >  sin*^  x2  ”  xi  *  ° 

Final  angle:  since  R  <  o  and 

*2  '  Z1 

z,  -  2^  >  o,  then - — - >  o 

and  P  •  +1 


Arc  length:  sin  *  c.p  <  sin  *  c_p  P  ■  +1 

1  o  2ro 

such  that  s,  -  Sj  >  o 


NOTE:  If  the  vertex  depth  were  contained  in  this  zone,  it  would  have  been 

used  as  the  lower  boundary. 
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The  progress  dlf ferentlates  between  each  of  these  four  cases  based 
on  the  gradient  and  initial  angle  at  the  starting  point.  The  appro¬ 
priate  boundarv  Is  then  established  (vertex  point  or  zone  boundary)  and 
the  ray  la  traced  to  this  boundary.  All  ray  parameters  are  evaluated 
and  a  test  for  convergence  Is  made.  If  the  ray  has  not  yet  converged, 
the  ray  trace  loop  Is  repeated  for  the  next  zone. 

When  the  ray  has  converged  on  the  surface,  bottom,  or  time  specified, 
all  of  the  parameters  necessary  for  evaluation  of  equation  (18)  In 
Chapter  II  are  saved  for  entry  Into  subroutine  PLOSS. 


2 •  Propagation  Loss 

The  propagation  loss  for  a  rav  traversing  N  zone  boundaries  is 

given  by 

2  2* 
r  cos  « 

L - = - 2 - ,2  ....  - -  (12) 


!N  x  -  x  ,  (x-x.,)  sln«  j 

V  n  n-1  +  * _ o  I 

'o  *-*  sin*  sin*  ,  sin:  sin#  { 
,  n  n- 1  n  f 

n«l 


*  ln< 


where  x  and  :  are  the  horizontal  coordinate  and  final  angle  at  each 
n  n 

zone  boundary  as  generated  by  RAYTRK;  subroutine  PLOSS  Is  used  to  evaluate 

equation  (12)  at  each  zone  boundary  and  at  each  convergence  point.  PLOSS 

also  takes  appropriate  action  for  the  following  special  cases. 

a.  If  the  ray  starts  with  initial  angle  4  '0,  then  equation  (12) 

o 

Indicates  that  all  zones, except  the  zones  In  which  the  ray  vertexes,  will  be 
characterized  by  the  following  propagation  loss. 

L  -  - 
x 
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In  a  vertex  rone  equation  (12)  la  valid  until  the  ray  actually  reaches 
the  vertex  point  at  which  point  the  Intensity  Is  Indeed  Infinite.  The 
subroutine  tests  the  final  angle  PHI2  for  zero  and.  If  true,  the  loss 
Is  set  equal  to  zero,  and  an  Indicator  Is  set  equal  to  Indicate  the 
infinity  condition.  Following  the  vertex,  equation  (13)  Is  used  until 
the  next  vertex  point  Is  reached. 


b.  For  vertexing  ravs  other  than  the  ray  with  zero  Initial 
angle,  we  use  the  limiting  form  for  the  propagation  loss  as  *  approaches 
zero  which  gives 


2  2« 
r  cos  t 
o _ o 

|(x-x^.)  sln«o 
•  ln*s 


(1*) 


The  program  makes  the  assumption  that  a  rav  will  never  be  allowed  to 
have  a  simultaneous  vertex  and  zone  transition,  l.e..  If  the  vertex 
depth  and  the  zone  depth  are  identical,  then  the  ray  is  assumed  to  vertex 
in  that  zone. 

At  the  conclusion  of  PLOSS  the  spreading  loss  has  been  computed 
for  one  of  rays  identified  by  TIMER.  Additional  subroutines  are  used 
to  supply  the  corresponding  surface,  bottom  and  attenuation  losses. 
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KAYTRK  Is  capable  of  selecting  from  a  set  of  four  absorption 
models  contained  In  subroutine  AL1.  The  two  models  presently  available 
are  zero  absorption  and  the  Thorpe^1'  model  where 


ATLOSS 


2  2 
O.lf  +  40 f 

1+f2  4100  ♦  f2 


db/kvd 


with  f  "  frequency  In  kilohertz 

4 .  Bottom  Loss 

RAYTRX  can  call  on  subroutine  BL1  for  either  zero  loss  or  an 
acoustic  province  and  Incident  angle  dependent  loss. 


5 •  Surf ace  Loss 

RAYTRX  has  provision  for  a  surface  loss  model  through  sub¬ 
routine  SL1.  At  present  this  choice  Is  limited  to  either  zero  or  3  db 
loss  per  bounce. 


Thorpe,  W.  H.f  Analytical  Description  of  the  Low  Frequency  Attenuation 
Coef f lc lent ,  J.  Acoustical  Soc  of  America,  Vol  42,  No.  1,  July  1967, 
p.  270. 
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CHAPTER  IV 


THE  TIMER  ALGORITHM 


In  Chapter*  II  and  III  the  mechanic*  of  ray  tracing  and  the  calcu¬ 
lation  of  propagation  loss  have  been  treated  In  some  detail.  This 
chapter  will  be  concerned  with  the  application  of  these  calculations 
through  the  use  of  the  TIMER  algorithm  to  the  evaluation  of  the  rever¬ 
beration  equation. 

The  reverberation  equation  contains  only  two  ray  parameters  explicitly, 
these  being  the  round  trip  propagation  time  and  the  propagation  loss  along 
the  round  trip  path.  However,  Implicit  in  this  equation  is  the  assumption 
of  spatial  continuity  of  the  round  trip  reverberation  path.  This  assump¬ 
tion  Is  of  some  practical  Importance  when  one  attempts  a  numerical  evalua¬ 
tion  of  the  reverberation  equation,  since  blstatic  reverberat Ion  paths 
( 1 . e .  ,  paths  for  which  the  outbound  and  inbound  ray  paths  are  different) 
are, In  manv  cases,  quite  significant.  Thus  we  seek  a  method  of  character¬ 
ising  each  ocean  environment  to  be  Investigated  in  a  manner  which  will  be 
consistent  with  the  assumpt ions  of  spatial  and  temporal  continuity.  The 
TIMER  algorithm  meets  these  requirements,  and  Its  operation  is  described 
as  follows. 


1.  First,  consider  the  problem  of  evaluating  the  propagation 
time  along  a  given  ray  path.  In  differential  form  we  have 


dt  1 

ds  ’  c(i) 


(1) 
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where  t  is  time,  s  Is  arc  length,  and  c(z)  Is  the  propagation  velocity 
along  the  ray  path.  Equation  (1)  Is  placed  In  a  convenient  form  for 
Integration  by  the  following  manipulations: 


dr _ 

clz)  sin,* 


c(z)  /l-cos' ♦ 


c(z)  /l-c(z)*p- 


where 


,  ?  •  — t  .  a  constant  for  a  given  ray. 

c(z)  c(z  )  ' 

o 


for  the  particular  case  in  which  c(z)  is  a  linear  function  in 
z,  equation  (2)  may  be  integrated  in  closed  form,  since 


c(z)  •  co  *  gz 


dz 


g  OR  dz  •  —  dc(z) 

—  ft 


f  rom  which 


t  .  t  .if**  _ tSlil _ 

1  8  •  c (z)  e'l-p-V 


1  c2  1  +Vl-C12p2 

- log  - 

*  e  ci  i  ♦*/; — n 

Vl-c2  p 


where  the  variation  in  z  from  to  z,  is  that  encountered  in  traversing 
the  ray  path  specified  by  p  from  Sj  to  s,.  Implicit  in  this  expression 
is  the  condition  that  the  integration  path  may  not  pass  through  a  vertex 
point,  since  the  integrand  is  unbounded  at  such  a  point.  However,  the 
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integral  is  convergent  in  the  limit  as  the  arc  length  approaches  a  vertex 
and  hence,  the  Integral  may  be  evaluated  in  steps  between  each  vertex. 

In  addition,  the  sign  to  be  used  with  the  radical  is  determined  by  noting 
that 

j: — T2 

\  1  -  c  p  •  sin* 
where  t  is  the  ray  angle. 

2.  The  basic  method  of  applying  equation  (3)  to  the  character¬ 
isation  of  a  given  ocean  environment  starts  by  tracing  from  the  source 
some  large  number  of  rays  distributed  in  equal  Increments  of  initial 
angle.  During  the  course  of  tracing  each  ray,  the  program  stores  the 
ray  parameters  at  a  set  of  predetermined  times.  These  times  usually 
correspond  to  the  times  at  which  reverberation  spectra  are  to  be  evaluated 
The  convergence  of  each  of  the  rays  to  the  desired  propagation  time  is 
accomplished  by  use  of  a  Newton's  method  subroutine  which  provides 
numerical  solution  of  equation  (3)  in  terms  of  the  rav  depth  z  for  each 
ray . 

The  result  of  this  effort  is  a  series  of  TIMER  tables  each  correspond 
ing  to  a  given  time  and  each  containing  a  complete  set  of  rav  parameters 
for  that  propagation  time.  Thus,  the  individual  TIMER  tables  constitute 
a  description  of  the  wave  front  at  each  specified  time  in  terms  of  the 
initial  ray  angles.  This  feature  of  the  TIMER  tables  allows  their  direct 
use  in  evaluating  the  reverberation  equat ion, s lnce  the  initial  ray  angle 
is  in  fact  sufficient  to  completely  describe  the  ray  path. 
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To  further  Illustrate  this  point  consider  the  graph  shown  in 
Fig.  IV-i .  This  particular  graph  is  a  plot  of  the  TIMER  data  correspond¬ 
ing  to  a  propagation  time  of  seven  seconds  for  a  bundle  of  rays  distributed 
in  initial  angle  from  -5  degrees  to  10  degrees.  The  source  is  located  at 
a  depth  of  10  yards  in  an  ocean  characterized  by  the  velocity  profile 
shown  in  Fig.  I V-2.  The  data  is  plotted  as  ray  depth  vs  initial  angle 
where  positive  initial  angles  represent  rays  with  a  downward  initial 
direction.  The  small  numbers  in  brackets  represent  the  number  of  surface 
and  bottom  bounces  taken  by  the  rays  which  comprise  each  branch  of  the 
TIMER  plot. 

As  an  example  of  the  application  of  this  data,  consider  a  volume 
scattering  layer  located  at  a  depth  of  100  yards.  Also  assume  that  a 
round  trip  propagat ion  time  of  fourteen  seconds  is  desired.  We  note  that 
the  line  drawn  across  the  TIMER  plot  at  a  depth  of  100  yards  intersects 
the  data  at  five  locations.  By  Interpolating  between  initial  angles,  a 
reasonably  good  estimate  may  be  made  of  the  initial  angles  required  to 
converge  the  rays  to  the  100  yard  depth  with  propagation  time  seven 
seconds.  A  ray  plot  of  these  five  rays  is  shown  in  Fig.  IV-3  along  with 
their  respective  propagation  losses  in  Fig.  IV-A.  This  method  is  suffi¬ 
cient  for  the  evaluation  of  the  reverberation  equation  when  both  the  out¬ 
ward  and  return  ray  paths  are  identical.  However,  when  a  bistatic  situation 
exists,  a  somewhat  more  involved  method  of  evaluating  initial  angles  must 
he  used. 

The  method  described  above  may  be  extended  to  the  case  of  a  plane 
surface  at  a  fixed  depth  cutting  through  a  number  of  TIMER  plots  as  shown 
In  Fig.  IV-5.  If  the  TIMER  tables  have  been  evaluated  at  close  enough 
intervals,  typically  one-half  second,  it  will  be  possible  to  generate 
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lines  of  initial  angles  vs  propagation  time  holding  the  depth  fixed. 

The  fixed  depth  would  of  course  be  the  surface  or  bottom  or  the 
boundaries  of  a  volume  scattering  layer.  This  data  forms  half  of  the 
solution  to  the  blstatlc  or  multipath  reverberation  problem.  However, 
we  still  need  a  method  by  which  to  relate  the  propagation  range  to  propa¬ 
gation  time,  since  the  blstatlc  path  must  be  coincident  in  range  and 
depth  as  well  as  having  the  same  round  trip  propagation  time. 

A  plot  is  made  of  range  versus  time  for  the  points  of  Intersection 
between  a  constant  depth  plane  and  a  sequence  of  TIMER  plots.  In 
Figure  IV-6  the  rays  chosen  are  the  direct  rays  1  and  2  of  Figure  IV— 1 
and  the  multibounce  ray  5.  Assuming  that  fourteen  seconds  after  trans¬ 
mission  is  a  time  of  interest,  we  may  reasonably  ask  if  there  is  any 
combination  of  dissimilar  propagation  paths  which  can  return  energy  to 
the  receiver  with  a  total  round  trip  propagation  time  of  fourteen  seconds. 

From  Figure  IV-6,  it  appears  that  there  is  indeed  such  a  combination 
of  paths.  The  requirement  that  both  paths  terminate  at  the  same  depth  is 
satisfied  by  the  cutting  of  the  TIMER  plots  at  a  fixed  depth.  The 
requirement  that  both  paths  terminate  at  the  same  range  is  satisfied  by 
constraining  the  range  variable  as  shown  in  Figure  IV-b.  Now,  by 
adjusting  this  range  variable,  R  ,  it  is  possible  to  find  a  combination 
of  propagation  times  tj  and  t,  whose  sum  is  the  required  fourteen  seconds. 
Knowing  the  propagation  times  tj  and  t^,  the  final  step  is  to  use  the 
graphs  generated  by  the  process  shown  in  Figure  IV-5  to  obtain 
Interpolated  values  for  the  initial  angles  of  the  rays  which  will  fulfill 
the  conditions  of  spatial  and  temporal  coincidence.  The  result  of 
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this  calculation  would  appear  similar  to  a  combination  of  rays  one  or 
two  and  ray  five,  as  shown  In  Figure  IV-3. 


The  process  of  multipath  identification  and  selection  has  not  as 
yet  been  Integrated  into  the  reverberation  program.  The  programs  for 
doing  this  task  are  under  development. 
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FIGURE  IV- 


FIGURE  IV  1 
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CHAPTER  V 


ARRAYS  AND  DOPPLER  BROADENING  CALCULATIONS 

I .  Array  Calculat Ions 

The  method  used  In  the  reverberation  program  for  the  evaluation 
of  the  hydrophone  arrav  directivity  is  basically  a  Fraunhoffer  diffraction 
model.  To  obtain  the  particular  form  of  the  directivity  equation,  we 
start  by  writing  the  pressure  resulting  from  a  single  transducer  as 

Pjfij.  t)  -  A ( J 5 1 1 )  «*p  J  ~  (jxj  -  ct)  (1) 

where  is  a  vector  from  the  1th  transducer  to  the  observation  point. 

The  function  A(  x^  )  is  an  amplitude  factor,  i  is  the  wavelength,  c  the 
propagation  velocity,  and  t  the  time.  The  total  pressure  due  to  an 
array  of  N  hydrophones  is 

N  , 

P< t )  -  £  A<  x  )  exp  J  —  (  x  -  ct)  (2) 

l-l 

At  this  point  we  Introduce  the  critical  assumption  that  the  distance 
from  a  reference  point  in  the  hydrophone  array,  sav  the  l1*’  transducer  to 
the  observation  point,  is  much  larger  than  the  distance  between  any  two 
elements  in  the  arrav.  With  this  assumption  we  mav  write  that  the  ampli¬ 
tude  functions  for  any  two  elements  are  very  nearly  equal,  since  their 
distances  to  the  observation  point  are  verv  nearly  equal.  Thus 

A(|b1|)  *  AflSjl)  (3) 
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and  extracting  the  constant  factor  exp  J  —  (-ct)  leaves, 

,  N  , 

P(t)  -  A(  x  )  exp(-j  y-  ct)  £  exp  J  ^  *, I 

1-1  1 


where  x^  Is  an  average  vector  from  the  array  to  the  point  of  observation. 

Let  the  Individual  vectors  x^  be  expressed  as  the  sun  of  the 

average  vector  x  and  a  location  vector  r.  for  each  element, 
o  1 


x.  ■  x  ♦  r. 
1  o  1 


Also,  let  ;  be  a  unit  vector  In  the  direction  of  observation,  l.e., 
parallel  to  the  vector.  Now,  using  the  original  assumption  that 
x^  Is  much  larger  than  ^  ,  we  can  write  the  following  approximation 
for  the  magnitude  of  x^  |  : 


V  *0  +  rl  '  ‘ 


Using  the  result  In  equation  (4)  gives 


P(t)  •  A( ! x  | )  exp  J  —  (l*0  -  ct)  exp  J  —  (r  •  *) 

1-1  ‘  1 

which  we  may  now  place  In  a  more  convenient  form  bv  normalizing  P(t)  In 
the  following  manner.  Define  P^^  as  the  maximum  of  equation  (4) 
obtained  by  holding  xq  constant  and  varying  <  over  4-  steradians. 

Thus  in  normalized  form 


£**P  J  —  <r.  •  t) 

Q Ct)  -  ^ - * - 


where  Q  is  the  maximum  of  the  summation  in  the  numerator. 
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It  Is  entirely  possible  to  associate  any  arbitrary  phase  with  each 


of  the  N  elements.  In  particular,  suppose  that  a  phase  equal  to 


2* 

'  IT 


IB  -  ♦  . 


(9) 


Is  associated  with  each  radiator  where  m  Is  a  unit  vector  in  some 


particular  direction.  Inserting  this  phase  In  equation  (8)  gives 


Q(4) 


exp  J  7~  (fj  *  ( t-m) ) 

1-1 


(10) 


Clearlv,  when  ■  >  m,  all  the  elements  are  in  phase  and  the  summation 
is  maximized  such  that  Q  -  N.  Equation  (10)  provides  a  computationally 

2)4  X 

convenient  method  for  accomplishing  plane  wave  phasing. 

In  addition  to  an  arbitrary  phasing,  we  may  also  allow  an  arbitrary 
weighting  function  to  be  applied  to  each  element.  For  most  practical 
purposes  a  weighting  function  dependent  on  the  angle  between  the  obser¬ 
vation  vector  and  a  unit  vector  characteristic  of  the  element  is 
sufficiently  general.  In  manv  cases  the  unit  vector  normal  to  the 
radiating  surface  is  used  as  a  reference  for  the  weighting  function. 

Taking  into  account  all  of  the  above  consideration,  we  mav  write 
for  the  general  case 


D (?)  -,^1- 


jrwfiij  •  *)  exp  j 


[r  <?, 


(*-■))  ♦  i 


•]f 


max 


(11) 


V-  3 


Arthur  9Aittlr)iu. 


where  Che  square  of  Q(i)  is  used  since  Che  direccivlcv  Is  normally 
associated  wlch  power  measurements. 

To  summarize  Che  variables  In  equation  (11)  are: 

th 

tv  unit  vector  characcerlsdc  of  Che  1  transducer. 

t  unit  vector  In  the  direction  of  observation. 


U(n.  •  t)  weighting  function  for  the  l*1’  element, 
v  location  vector  for  the  lll>  element. 


plane  wave  phasing  unit  vector. 


0 

max 


arbitrary  phasing  for  the  I1*'  element 

normalizing  value  equal  to  the  maximum  of  the 

*♦ 

numerator  over  all  values  of  i . 


As  an  example  of  the  use  of  equation  (11),  we  show  in  Fig.  1  and  Fig.  2 
the  vertical  and  conical  patterns  associated  with  an  array  of  hydro¬ 
phones  having  cylindrical  symmetry. 

I 1 .  Doppler  Broadening  Calculations 

The  reverberation  equation  requires  the  evaluation  of  a  convolution 
integral  of  the  form 


>(f) 


F(<j)  A(f-o)  dr 


where  F(-)  is  a  function  characteristic  of  the  scatterer  motion  and 


V-4 


(12) 


Srthur  21  \ittlr  >nr 


A(f-o)  is  Che  spectrum  of  Che  transmitted  signal.  In  general  the  repeti¬ 
tive  evaluation  of  this  integral  by  conventional  techniques  is  not  feasible. 
However, by  making  use  of  the  Fast  Fourier  Transform  (FFT)  technique,  this 
calculation  may  be  performed  with  relative  ease.  In  fact,  the  use  of  the 
FFT  greatly  simplifies  a  number  of  the  computational  aspects  of  the  rever¬ 
beration  problem. 

Before  investigating  this  simplification,  we  will  define  the  FFT  and 
its  inverse  and  give  without  proof  some  elementary  but  useful  properties. 

The  Discrete  Fourier  Transform  (OFT)  of  a  sampled  signal  of  limited  band¬ 
width  is  defined^  as 


F 

n 


fk  CXP 


(J  jf-  nk) 


(13) 


and  the  inverse  is 


f 


k 


N 

N  Fn  exp  ! T  nk) 

n»l 


(14) 


Using  conventional  techniques  the  evaluation  of  either  equation  (13)  or 

2 

equation  (14)  will  require  N  complex  operations  (i.e.,  a  sine  or  cosine 
look  up  and  a  multiply).  The  FFT  derives  its  time  savings  by  capitalizing 
on  the  high  degree  of  svraraetry  of  the  conversion  matrix 

W(nk)  ■  exp  J  jj-  nk  (15) 

in  order  to  greatlv  reduce  the  number  of  complex  operations  required. 

^Tlie  definition  of  the  DFT  is  not  uniform  in  the  literature.  Some  authors 
use  A^/N  as  the  DFT  coefficients,  others  use  A^/tN,  still  others  use  a 

positive  exponent.  Cochran,  W.  T.,  et  al.  What  is  the  Fast  Fourier  Transform? 
Proc  of  the  IEEE,  Vol  55,  So.  10,  October  1967,  p.  1665. 
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Among  the  useful  properties  of  the  DFT,  we  have  that  the  time 


average  power  In  a  sample  signal  {  f  ^ }  is  given  by 


P 


-  £ 

n- 1 


F  F 
n  n 


(16) 


In  addition,  we  note  that  the  Individual  components  of  the  power 
density  spectrum  of  the  sampled  function  (f^)  ere 

P  -  F  F  *  (17) 

n  n  n 


which  Implies  that  the  autocorrelation  function  of  the  signal  (f.  } 
may  he  obtained  by  inverse  transformation  of  the  set  of  samples  ( Pn } . 
That  is 

♦k  •  £  pn  exp  (_J  ir  nk) 

n«l 


(18) 


Another  useful  property  Is  that  convolution  In  either  frequency  or  time 
corresponds  to  multiplication  in  the  inverse  domain.  Thus  Che  evalua¬ 
tion  of  the  convolution  Integral  would  proceed  if  done  bv  conventional 
methods  as 

=„  •  £  s  Vk  <19> 

K  ”  1 

However,  by  making  use  of  the  above  properties  and  transforming  (G^) 
and  {A.  ),  we  obtain 

*n  "  £  Vk  exp  (J  r  nk)  (20) 

k-1 

which  Is  again  evaluated  bv  means  of  the  FFT. 
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The  sequence  of  operations  as  used  by  the  program  is  as  follows. 

The  input  consists  of  a  time  function,  f(t),  representing  the  trans¬ 
mitted  signal  and  a  set  ot  numerical  values  representing  the  variances 
of  the  Doppler  spreading  from  both  the  surface  and  the  volume  scatters. 
The  input  signal  is  sampled  by  a  signal  generator  subroutine  to  give  set 
of  sample  values. 

f(t)  ♦  sampled  ■*  {f^} 

This  set  Uk)  is  then  transformed  after  which  the  magnitude  squared  of 
the  transform  is  obtained. 

{f.  )  -  FFT  -  ( F  > 
k  n 

F  •  magnitude  squared  •  I A  ) 
n  n 

Now  the  set  P  }  is  the  power  density  function  to  he  convolved  with  the 
n 

scatterer  function.  Next,  we  generate  the  transform  of  the  scatterer 

characteristic  function  g^}  which  we  assume  to  be  gausslan.  Following 

this  operation  the  spectra  \  )  is  transformed  to  give  the  sampled 

n 

correlation  function  >,  i. 

k 

{A  )  -»  Inverse  FFT  -•  (A.  ) 
n  k 

Multiplication  by  the  scattering  function  { g^ )  followed  by  a  FFT 

returns  the  set  { 5  )  representing  the  power  density  spectra  of  the 
n 

Doppler  broadened  input  signal. 
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{>k)  •  (g^  *  “ultiplication  ♦  (4^) 


F igures 

case  of 


Wk>  -  FFT  -  { *n) 


through  8  Illustrate  this  sequence  of  operations  for  the 
linear  f.m.  pulse  using  the  following  parameters: 


sampling  interval 

pulse  length 

repetition  rate 

rate  of  change  of 
frequency  with  time 

total  signal  bandwidth 

bandwidth  increment 


0.005  sec. 

0.50  sec. 

1.00  sec. 

100.0  cps/sec. 
100  cycles 
0.5  cycles 
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CORRELATION  FUNCTION  DISPLACEMENT  TIME  ^  SEC. 
CHARACTERISTIC  FUNCTION  FOR  QOppLER  BROADENING 


CORRELATION  FUNCTION  DISPLACEMENT  TIME  -  SEC. 
PRODUCT  OF  THE  CHACTERISTIC  AND  CORRELATION  FUNCTIONS 


FREQUENCY  IN  CYCLES  PER  SECONO 
POWER  DENSITY  SPECTRUM  OF  THE  DOpPLER  BROROENED  SIGNRl 


CHAPTER  VI 


REVERBERATION  SPECTRUM  COMPUTATION 

The  reverberation  sped  run  time  historv  ia  generated  by  the  program 
RAND  which  almply  acta  aa  a  driver  for  all  of  the  auhroutlnea  described 
in  the  preceding  chapters.  A  series  of  times  is  chosen  at  which  spectra 
are  desired.  Programs  TIMER  and  TIMERP  are  run  for  these  times  with  the 
specified  ocean  and  source  depth.  Prom  the  TIMER  plots  or  printed  output 
initial  angles  are  obtained  for  each  reverberation  path  (first  order 
surface,  first  order  bottom,  etc.)  for  each  time.  This  step  also  allows 
an  estimate  of  the  intensity  of  the  selected  ravs ,  since  the  space 
between  points  on  the  TIMER  plot  is  proportional  to  the  spreading  loss. 

Program  RAPID  is  then  run  to  produce  plots  of  the  array  directivity. 
Since  RAND  has  to  integrate  equation  1-40  over  all  azimuth  and  depres¬ 
sion  angles,  computation  time  can  be  decreased  if  it  is  not  necessary 
to  Integrate  over  all  of  the  side  lobes  in  detail.  Using  the  RAPID 
plots  one  can  select  integration  limits  within  which  the  array  pattern 
is  computed  in  detail  and  outside  of  which  a  constant  sidelohe  level  is 
specif  led. 

Doppler  broadening  characteristic  functions  are  specified  for 
program  RAND  by  estimating  the  bandwidth  of  the  scatterer  motion  spectrum. 
The  ocean  is  specified  by  its  velocity  profile,  the  surface  bv  the  wind 
velocity  and  the  bottom  by  its  acoustic  province.  Two  types  of  trans¬ 
mitted  signal  are  presently  available,  a  pulsed  CW  signal  and  a  linear 
FM  signal. 


c 

0 
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Figures  6.2  and  6.3  show  the  output  of  the  program  for  the  case 


of  a  0.123  second  duration  CW  pulse  whose  power  density  spectrum  Is  shown 
in  Figure  6.1.  The  particular  reverberation  paths  used  in  this  example 
were  determined  by  examination  of  the  TIMER  program  output  for  the  one 
way  propagation  times  of  14.0  and  13.0  seconds.  Within  the  depression 
angle  range  considered,  it  was  found  that  three  rav  paths  existed  for 
the  14  second  case  and  four  paths  existed  for  the  13  second  case.  The 
first  three  paths  consisted  of:  1)  direct  to  the  bottom,  2)  one  bottom 
bounce  and  terminating  In  a  scattering  layer,  and  3)  one  bottom  bounce 
and  terminating  on  the  surface.  The  13  second  time  allowed  one  addi¬ 
tional  path  consisting  of  one  bottom  bounce,  one  surface  bounce  and 
termination  in  the  scattering  layer.  The  initial  angles  of  each  of 
these  rays  now  formed  the  basic  Input  for  each  reverberation  time. 

The  transducer  array  used  In  this  example  was  depressed  an  angle 
of  30*  and  the  following  doppler  broadenlngs  were  assumed  for  the 
various  scatterers: 

ship  1  cps  rms 

surface  2  cps  rms 

volume  4  cps  rms 

For  this  particular  example  it  happens  that  the  first  path,  direct  to 
the  bottom  and  back.  Is  dominant.  For  other  times,  oceans,  nrravs  and 
depression  angles,  the  dominant  reverberation  may  come  from  a  less 
obvious  path.  The  RAND  program  allows  a  svstematlc  investigation  of 
the  effect  on  reverberation  of  sonar  and  ship  operating  parameters. 
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FREQUENCY  IN  CYCLES  PER  SECONO 

TRANSMITTED  SIGNAL  SPECTRUM 
REVERBERATION  TEST  USING  VP-H  WITH  SOURCE  AT  100  TAROS 

REVERBERATION  TIRE  1.00  JUNE  21.1969 
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FIGURE  VI-2  SPECTRA  FROM  INDIVIDUAL  REVERBERATION  PATHS 
AND  TOTAL  POWER  SPECTRA  FOR  ROUND  TRIP 
PROPAGATION  TIME  OF  28  SECONDS 
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FIGURE  VI-3 
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SPECTRA  FROM  INDIVIDUAL  REVERBERATION  PATHS 
AND  TOTAL  POWER  SPECTRA  FOR  ROUND  TRIP 
PROPAGATION  TIME  OF  30  SECONDS 


